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Abstract. 

We compare multiplicities as well as rapidity and transverse momentum distributions of 
protons, pions and kaons calculated within presently available transport approaches for heavy 
ion collisions around 1 AGeV. For this purpose, three reactions have been selected: Au+Au at 
1 and 1.48 AGeV and Ni+Ni at 1.93 AGeV. 



1. Introduction 

Heavy ion collisions in the energy range from 50 AMeV to several AGeV are quite complex. The 
single particle spectra are not isotropic and the anisotropy depends significantly on the ccntrality 
of the reaction. With increasing energy one observes the rise and fall of multifragmcntation, 
where multifragmentation means that in each single collision a multitude of intermediate mass 
fragments is produced. While above Ekin = 289(1583) A/ey, 7r's(A'"*"'s) can be produced in 
elementary nuclcon-nucleon collisions, experiments have shown that in collisions between nuclei 
Fermi motion and multiparticle effects allow for meson production at energies well below the 
threshold in an elementary reaction. Finally nucleons can be excited to resonances which - 
according to observations in 7 A reactions - change their properties in a nuclear environment. 

This complexity of the reactions makes it very difficult to link the experimental results in 
a unique way with the specific underlying physical processes. The challenge to interpret the 
experimental results was first met by Yariv and Frankel and Cugnon [5] who took advantage 
of the increasing computer power at that time and developed the so-called "cascade approach" 
in which nucleons behave similarly to classical billiard balls and scatter with the free elementary 
cross sections. This computational approach allowed for completely new insights into the reaction 
mechanism and allowed for the first time to interpret the nonequilibrium features of the observed 
spectra. Disregarding binding energy and mean field effects, these approaches have been limited to 
small systems at higher energies. The extension to lower energies became possible half a decade later 
with the numerical realization of the Boltzmann-Uchling- Uhlcnbeck (BUU) ^ or Vlasov-Uehling- 
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Uhlenbeck (VUU) 0j approaches which supplemented cascade approaches by the introduction of 
an attractive mean field and of Pauli blocking, which suppresses collisions if the phase space is 
already occupied by other nucleons. This made it possible to describe heavy ion reactions down 
to energies of 20 AMeV. Over the years these approaches became more and more refined as a 
result of constructive competitions among different groups and many physical questions have been 
addressed and answered. The pre-equilibrium proton emission has been quantitatively reproduced. 
It is nowadays established that subthreshold mesons like r], K or p are created in elementary 
collisions, where at least one of the collision partners had gained additional energy in previous 
collisions. It is further established that resonances are the primary source for the produced pions. 

In the meantime, the simulation programs have matured to the level which allows to go 
even beyond the original goal of describing the elementary reaction features: They are nowadays 
commonly used as a tool to find out whether the experimental results allow for a determination of 
physical quantities for which no solid theoretical predictions are available. These quantities include 
the energy which is needed to compress nuclear matter as well as the properties of resonances and 
mesons in a hadronic environment. 

For the following three reasons it seems to be useful to compare the simulation programs which 
have been developed over the years: 

1) The simulation programs are rather complex - consisting of several thousand program lines 
- and use quite different numerical techniques. Thus, one has to assess whether the different 
numerical procedures lead to the same results. 

2) All programs rely on inputs. Those inputs include all the needed elementary cross sections. 
Among those there are cross sections which are not known experimentally (like those including a 
baryonic resonance in the entrance channel), and different theoretical implementations have been 
introduced. It is evident that different elementary cross sections can yield different results. 

3) In all programs resonances are produced. Their properties in the hadronic environment as 
well as their time evolution are little known. The programs use different parameterizations for this 
and it is not a priori evident how those differences influence the final result. 

Thus the detailed comparison presented in the present paper may serve as the first step 
towards a critical assessment of the predictive power of the simulation programs and to identify 
the features which have to be improved. We would like to mention that some elements towards 
a future development are already available. Very detailed calculations of the spectral function 
of kaons in infinite matter have been published [S] and first steps towards describing the time 
evolution of resonances in matter have been made |S1|S1[7|- These advances, based on the gradient 
expansion of the Kadanoff-Baym equation, are difficult to formulate for test particles, and further 
approximations have been shown to be necessary for actual applications 8 . Apart from these 
approaches there are a couple of studies providing information on the off-shell transition rates 

pun]. 

2. The simulation programs 

In this comparison nearly all the presently available simulation programs have taken part: The 
simulation program from the Budapest/Rossendorf group the HSD program developed in 
Giessen by Bratkovskaya and Gassing |12[ll3j . the simulation program developed by Bratkovskaya, 
Effenberger, Larionov and Mosel in Giessen ^J], the code developed by Danielewicz in Michigan 
|15| . the RVUU approaches of the Texas A&M group and the Munich/Gatania/Tiibingen group 
|17| . Besides these true BUU or VUU models, there are several models which are based on the 
Quantum Molecular Dynamics approach. These n-body approaches allow for the description of 
fragment production but as far as the observables discussed here they should give the same results. 
Therefore, we just refer to ref. JS] for the general features and do not discuss the differences 
here. These simulation models include the IQMD and URQMD [201, both developed in a 
collaboration between Nantes and Frankfurt, as well as the QMD of the Tiibingen group (by 
Fuchs) |21 l . These programs contain a different number of baryonic resonances as can be inferred 
from table The URQMD code has been developed for higher (SPS to RHIG) energies and 
includes neither a binding energy for nucleons nor a kaon-nucleon potential but all known mesonic 
and baryonic resonances. Also its production mechanism for is different from other approaches. 
Therefore it is only included to illuminate the differences. 
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Program 


r 


Resonances 


Barz 




A's,N's with M < 2 GeV 


RVUU (Chen) 


'^iso 


A(1232) 


Daniclcwicz 


Tiso 


A(1232),7V*(1440) 


QMD (Fuchs) 


'^iso 


A(1232),7V*(1440) 


Gaitanos 


'^iso 


A(1232),7V*(1440) 


HSD (Gassing) 


'^iso 


A(1232),7V*(1440),7V(1535) 


IQMD (Hartnack) 




A(1232) 


Larionov 


1/120 MeV 


A's, N's with M < 2 GeV 



Table 1. The width of the A resonance and the resonances employed in the programs. 

3. The Procedure 

The simulations programs, which are compared in this article, have been frequently compared with 
experimental data and usually gave quite satisfying agreement. It is, however, difficult to base 
program comparisons on those calculations because rarely the same data have been compared with 
different simulation programs and in the rare cases one has done so, it is not evident that the same 
experimental filter or the same centrality selection has been employed. To assess the predictive 
power of these programs and to see whether differences in the predictions are sufficiently large 
to be experimentally relevant, it is much better to compare the results of the programs directly 
without any cuts. In addition, all programs employ the impact parameter as an input variable and 
therefore it is most convenient to compare the results at a given impact parameter. This excludes 
of course any comparison with experiment. In this publication we concentrate on 3 reactions; Au 
+ Au at 1 and 1.48 AGeV and Ni+Ni at 1.93 AGeV, all at an impact parameter of 6 = 1 fm as well 
as the pion yield as a function of the impact parameter. This choice was guided by the available 
experimental results. 

4. Protons 

The final proton rapidity distribution measures the amount of stopping and reflects how much 
energy becomes available for particle production. The recent analysis of the FOPI collaboration 
j22| shows that, in the center of mass system, the variance of the transverse rapidity distribution 
is always smaller than twice the variance of the longitudinal rapidity distribution even in the most 
central collisions. If the system were in thermal equilibrium one would expect that the ratio R of 
the variances equals one. Experimentally the largest degree of thermalisation (R « 0.85) is observed 
for central Au + Au collisions around 400 AMeV. Below that energy Pauli blocking is still too 
active, while above the NN cross section becomes more and more forward peaked and therefore the 
energy transfer into the perpendicular directions is lowered. In addition, more energy has to be 
transfered in order to come to equilibrium. The rapidity and transverse momentum distributions 
from simulations, figs^andl^l reflect this approach to equilibrium. For the heavy systems the 
final rapidity distributions is peaked at or very close to midrapidity and the simulations differ 
only little. Initially the protons are located around ycm = 0.68,0.8,0.89 for the three reactions, 
respectively. The resulting rapidity shift for the large systems is thus large. A slightly smaller 
shift is observed for the light system and there also differences develop between the predictions of 
the different programs. This kind of differences should be measurable and can help to validate the 
assumptions made. In transverse direction and on logarithmic scale, we find a general agreement 
among the programs for Au+Au collisions. However, in the lighter system pronounced differences 
develop. They seem to be correlated with the rapidity distribution and reflect a different degree 
of stopping in the different approaches. The differences are most pronounced at small transverse 
momenta. To fill this part of phase space many collisions are necessary. Therefore its form is 
a measure for the number of collisions during the interaction. Table |21 summarizes the average 
squared momenta in the center of mass system at the end of the reaction, along the beam (z) 
direction and in the direction (y) perpendicular to the reaction plane. 
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Figure 1. Final proton rapidity distribution at b = Ifm in the different approaches 
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Figure 2. Final proton transverse momentum distribution at 6 = Ifm in the different approaches 





1 AGeV Au+Au 


1.48 AGeV Au+Au 


1.93 AGeV Ni+Ni 


Program 


{pD 




Tpf) 


{pD 


^pD 


W) 


ipl) 








(GeV'/c') 




[GeV'/c') 




{GeV'/&) 




Barz/Wolf 


0.157 


0.109 


0.69 


0.230 


0.146 


0.64 


0.351 


0.147 


0.42 


RVUU (Chen) 


0.165 


0.121 


0.73 


0.256 


0.172 


0.67 


0.416 


0.172 


0.42 


Danielewicz 


0.162 


0.102 


0.63 


0.257 


0.126 


0.49 


0.407 


0.104 


0.26 


QMD (Fuchs) 


0.158 


0.111 


0.70 


0.262 


0.138 


0.53 


0.413 


0.107 


0.26 


Gaitanos 


0.131 


0.116 


0.89 


0.240 


0.142 


0.60 


0.300 


0.104 


0.35 


HSD (Gassing) 


0.138 


0.120 


0.87 


0.219 


0.157 


0.72 


0.361 


0.148 


0.41 


IQMD (Hartnack) 


0.143 


0.121 


0.85 


0.226 


0.161 


0.72 


0.378 


0.158 


0.42 


Larionov 


0.163 


0.124 


0.76 


0.252 


0.159 


0.63 


0.405 


0.136 


0.34 



Table 2. Average out-of-plane and longitudinal squared momenta for protons in the different 
reactions. 



5. Pions 



The direct pion production cross section in nuclear interactions is small. Therefore all pions are 
created in simulations via resonances. These resonances are produced with the free inelastic NN 
collisions - which is well known - and in the region of interest comparable to the elastic cross section. 
In our energy region these are mostly A resonances which disintegrate into a tt and a nucleon. The 
problem is how to treat this resonance in matter as is required for the simulation programs. The 
A resonance, as noted in the particle data booklet, has a width F of about 120 MeV and hence 
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Figure 3. The different A lifetimes wliicli have been proposed. 



a lifetime of t = i = 1.7/m/c. Therefore it disintegrates inside the nuclear medium. The width 
has been determined by analyzing the phase shift in ttN scattering as a function of ^/s. Around 
S{^/s) = 7r/2 the dependence on the energy may be approximated with S{^/s) = tan^^{ 2i^ml- ^/i) ) • 
The lifetime as inverse of a constant width V is justified if a) the width were small as compared to 
the mass difference between the resonance and its decay products (F << mA — (toat + m^)) and if 
b) the energy spread of the particle wave packets were large compared to F. This is, however, not 
the case: the nucleons are represented by delta functions in energy. Therefore one has to address 
the question which lifetime should be assigned to the object which is created in an inelastic nucleon- 
nucleon collision with a sharp ^/s and which yields two nucleons and a pion in the exit channel. 
There are two approaches in the literature. Ericson and Weise in particular, have analyzed 
the ttTVA system in a relativistic isobar model and found for the scattering amplitude 



/33 - _2 „ • .11,1,1 ' (1) 



TO^ — s — i7|k|'^ 



a form which resembles the usual relativistic Breit-Wigner form and yields an energy dependent A 
decay width of 

With /a = 4.02 this form reproduces very well the tt N data. rA(|k|) decreases with decreasing 
^/s and therefore a lifetime defined as 

(3) 



FA(|k|) 

increases with decreasing y^s. 

On the other hand, the lifetime of a resonant state has been already analyzed by Wigner . 
In recent time this considerations have been generalized by Danielewicz and Pratt |25j . They found 
that quantum mechanically the reaction slows down with a time delay between the outgoing and 
the incoming wave of 

d5{E) 



».s-2 . (4) 

Hence maximal time delay occurs at the center of the resonance. Fig. |31show the A lifetimes in the 
different approaches as a function of the A mass. Randrup, Kitazoe and Huber are synonymous for 
three rather similar approaches which calculate the width using the p-wave scattering amplitude 
as described above. The description used in the standard versions of the different programs are 
shown in Tabled 
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The influence of the choice of r on the final number of pions is by far not neghgible. Almost 
independent of energy and size of the system Tiso may yield up to 60% more pions than r^^g . If the 
A has a short lifetime, several generations of tt and A are produced. In an expanding system the 
average available center of mass energy in the NN system decreases and the cross section for the 
production of new A's becomes smaller whereas for the backward reaction it is little influenced. 

Fig. 01 shows the number of produced pions as a function of the impact parameter for the 
three reactions. We see the expected significant variations of the total yield. In order to put the 




4 8 12 
b [fm] 



4 8 12 
b [fm] 



Figure 4. Final 7r+ and vr multiplicity as functions of the impact parameter in the standard 
versions of the different approaches characterized by different descriptions of the A lifetime (see 
table EJ. 



results on an equal footing we have decided that all other calculations presented here use a constant 
lifetime of r = l/r,r = 120 MeV. The rapidity and transverse momentum distributions, obtained 
with this constant hfetime, are displayed in figs. |5landEl respectively. We see, first of all, that 
the differences between the calculations are reduced but not eliminated. The calculation of Reiter, 
based on URQMD, has to be interpreted with care because in this calculation neither potentials 
are present nor - being an approach for high energies - special care has been taken to model the 
inelastic cross section close to threshold. In the program of Chen only the isospin averaged pion 
yield is available. For that program we have therefore plotted here for tt~ and tt"*" one third of 
the total TT yield. Aside from Reiter's results, the rapidity distributions are rather similar for the 
majority of the programs. There are, however, some exceptions which show a narrower distribution. 
The experimental ratio R — 7r~/7r+ at 1 AGeV is close to the value R=1.95 which one expects 
from the Clebsch Gordon coefficients if (in this neutron/proton asymmetric system) the tt's are 
exclusively produced via A resonances. Most of the simulation programs reproduce this value. The 
slopes of the transverse momentum spectra at midrapidity are very similar as well. The differences 
concern exclusively tt's with very small pt values. This part of the spectra is sensitive to A's 
with a very small mass. We would like to stress that some of the programs include higher {N* or 
A*) resonances whereas others are limited to the A resonance. This may be of importance at the 
higher beam energies. Generally the simulation programs seem to overpredict the pion yield. A 
compilation of the available data in ref. | 26| shows an average total pion yield per participant of 
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Figure 5. Final 7r+ and n rapidity distributions at 6 = Ifm in the different approaches with 
an enforced A Hfetime of 1/120 MeV. 




Figure 6. Final tt"*" transverse momentum distribution at b= 1 fm and |j/cm < 0.5| in the 
different approaches with an enforced A lifetime of 1/120 MeV. 



TT^ + TT^ + TT" = 1.5(7r"'" + 7r~) = 0.075 and therefore about 30 pfons in the most central colHsions of 
Au+Au at 1 AGeV. Please note that the slope of the pion transverse momentum spectra is rather 
different from that of the protons. This is due the convolution of the A transverse momentum 
distribution with its the decay kinematics. 



6. K+ Mesons 



Fortunately, the uncertainty about the optimal description of the A lifetime influences neither the 
multiplicity nor the form of the spectra in a significant way because low mass A's have a 
small cross section to produce a K^. Due to the high threshold in the mcsonic channel the 



Transport Theories for Heavy Ion Collisions in the 1 AGeV Regime 



8 



are mostly produced in baryon-baryon collisions. To make comparisons more straightforward we 
enforced for the comparison the same cross sections for the NA — > and the AA 
channels in the different codes. These cross sections are dominant at Ekin > 1.5AGeV. Therefore, 
and due to the constant A lifetime, the results presented here may differ from already published 
results. For the NN cross section we kept the original cross sections of the programs. 

For the channel pp new data at low energies are available [2HI- In addition there is the 

problem of how to extrapolate this pp cross section to the nn, np — > if + channels. Usually this 
extrapolation is done by assuming that the K'^ is accompanied by a A and that only the isospin 
coefficients change. In the case of the this is, however, not so easy because this extrapolation 
depends on whether a kaon or a pion is exchanged. For the case of a pion exchange we find 
a{np K^) = 5/2 <T{pp K^) and consequently a{NN K^) = 3/2 a{pp K^); in the case 
of a kaon we find (T{pp — > K~^) = 2a{np and therefore a{NN K^) — 1/2 (j{pp K~^). 

Very recent data [211 point, however, towards a a{np K^)/a{pp —* K^) ratio of 3-4. A further 
complication is the production of the in the pp K^NY, channel, which has been found 
to be small |2H1j but this might be a consequence of the final state interaction. The different 
a{BB K~^) cross sections employed are presented in fig. [7| and Table |2| Top panels of 



Program 


a{NN K+) 


cr(iVA - 


K+) 


(aAA -> K+) 


Barz 


M 


1^ 







Chen 


27| 


27 






27 


Fuchs 


30^ (A) S 


27 






27 


Gaitanos 


30 (A) 27 S 


27 






27 


HSD (Gassing) 


.12i 


27 






27 


IQMD (Hartnack) 


30 (A) EH S 


27 






27 


Larionov 


123 


|27, 






27 



Table 3. The K'^ cross sections used in the programs. 




Figure 7. The isospin averaged cross sections NN — > . 

fig. and fig. |3 present, respectively, the rapidity and transverse momentum distribution 
obtained for those cross sections. Also here the URQMD calculations (Reiter) cannot directly be 
compared with the others. In this approach, the are created by two body resonance decays 
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Au+Au 0.96 AGeV 



no potential 




Figure 8. Final rapidity distribution at 6 = 1 fm and with an enforced A lifetime of 1/120 
MeV (top row without, bottom row with KN potential) in the different approaches. 



whereas in all other approaches the production is a process with three bodies in the final state: 
BB K+A{J:)B. We see that URQMD is not suited for 1 AGeV but aheady at 1.5 AGeV it 
describes the spectra in a reasonable way although spectra measured in elementary collisions point 
towards a three - body process. In order to understand the 50% differences of the predictions 
one has to realize that the are produced at subthreshold energies in the Au+Au collisions. 
Therefore the nucleons, which are involved in the production, had already collisions before. 
This encodes in the production yield the whole history of the reaction. In addition, the 
production cross section rises exponentially close to the threshold and therefore small differences 
in the parametrization of the Fermi motion of the nucleons or in the potential of the hyperon change 
the cross section considerably. Some programs contain, in addition, the channel ttB —^ YK~^, but 
usually this channel does not contribute a lot at these energies due to its large threshold and the 
few pions as compared to baryons. 

Kaons change their properties in the medium. This is the prediction of many theoretical 
approaches which range from chiral perturbation theory to Nambu Jona-Lasinio model 
calculations 31 . All these theories should give reliable predictions if the density is low but suffer 
from the fact that the extrapolations to higher densities are less and less well founded. For our 
purpose, it is sufficient to parametrize the density dependence by 

ujK{pB,k) = \J m*j^{pBY + (5) 
where the effective mass is 

= ( 1 - a— ) , (6) 
V PoJ 

with pb and po being the baryon density and the normal nuclear matter density, respectively. 
Table 0] displays the different values of a employed in the calculations: It has been suggested in 
ref. 32 to compare the ratio '^^^^ I '^dpt ' '^^li*^'^ been measured [HHIj with calculations using 
different K^N potentials 33 , to fix the potential at normal nuclear matter density. 
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Figure 9. Final transverse momentum distribution at 6=1 fm, \ycm, < 0.5| and witli an 
enforced A lifetime of 1/120 MeV (top row without, bottom row with KN potential) in the 
different approaches. 



Program 


a[MeV] for K+ (R-) 


Barz 


-0.05 (0.16) 


Chen 


-0.04 (0.22) 


Fuchs 


-0.07 


HSD (Gassing) 


-0.04 (0.10) 


IQMD (Hartnack) 


-0.075 (0.22) 


Larionov 


-0.06 (0.2) 



Table 4. Parameters for the density dependence of the in medium kaon masses m'^ ("ealGl. 



Despite of the moderate increase of the u!k+ {k = 0) with density, the potential changes the 
yield considerably, as one sees in figs. |Hland|5| This is due to the fact that kaons are produced 
at high density when the mean-free-path is short and therefore the A's have a higher chance to 
collide with a nucleon. Depending on the density dependence of the potential, we see a decrease 
of the total yield of about a factor of 2 for the heavy systems and less for the light system if 
the potential is switched on. The form of the transverse momentum distribution by Fuchs differs 
from those from other calculations because he uses a covariant form of nuclear current instead of 
only its forth component, the baryon density, as the other approaches. 

The elementary creation process is azimuthally isotropic. It has been argued that the 
azimuthal distribution of the contains information on the nucleon potential [34|. This is 
based on the fact that emitted into the azimuthal angle (j) = 90° come from highly compressed 
matter whereas those being emitted near — 0°(180°) have traversed the cold spectator matter. 
Therefore one expects an enhancement around (j) — 90° due to the more repulsive potential. Such 
an effect may, however, be superimposed by the consequence of the different path lengths in matter. 
Fig llOl shows the result. The majority of the programs produce a different azimuthal distribution 
without (top row) and with (bottom row) K^N potential. Without the potential, the distribution 
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is more isotropic and the expected enhancement is seen if the repulsive potential is active. 
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Figure 10. Final azimuthal distribution at b=5 fm, \ycm < 0.5| and with an enforced 

A lifetime of 1/120 MeV (top row without, bottom row with KN potential) in the different 
approaches. 



7. K- Mesons 

The production process for the K~ mesons in nuclear reactions is much more complex than that 
for the K^. The reason for this is their s quark content. On the one hand it is much harder to 
produce a K~ in elementary pp reactions because of its higher threshold {Exin — 2.513 GeV for 
pp K^K^pp), and on the other hand a K~ can exchange its s quark with strange baryons 
{K~N A(I])7r). The cross sections for these exothermic reactions are very large. The inverse 
reaction A(E)7r — > K~N has been identified as the dominant source of the K~ production in 
heavy ion reactions at the energies we consider here |35| . This process is absent in elementary 
pp interactions and therefore the K~ yield in heavy ion collisions is orders of magnitudes larger 
than that in pp collisions for energies not too far away from the threshold. This process couples 
the to the A(I]) and hence to the production [211 EZI ISH] and explains therefore the 
experimental observation "39' that the to K~ ratio is independent of the impact parameter 
The K~ production in baryon-baryon interactions suffer from the large threshold and takes place 
only early when the system is dense. Thereafter the K" gets easily reabsorbed. In figs. IlllandlT^we 
display the rapidity and transverse momentum distributions. The URQMD calculations suffer 
from the fact that they produce too few strange baryons (produced in the BB K^Y reactions) as 
can been seen from the low multiplicity. The channel confronts the simulation programs 
with the challenge to calculate the K~ N and the Att cross sections at values of ^/s which are in 
free space below the threshold. These problem has been solved differently by the different groups. 

In view of the complexity of the production processes, the results of those codes, which 
can be directly compared, are rather similar, especially at higher energies. At higher energies, 
the transverse momentum distributions are also rather similar. At lower energies the HSD code 
produces a somehow softer distribution. Also the K~ changes its properties in the medium. As 
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Figure 11. Final K rapidity distribution at b=l fm and witli an enforced A lifetime of 1/120 
MeV (top row without, bottom row with KN potential) in the different approaches. 



pointed in 001) due to the strange resonances this change is rather complex. In fact the detailed 
calculations [H] , based on the well constrained KN interaction from 01] , raise the question whether 
the can still be treated as a quasi particle at high densities and temperatures. On the other 
hand, most of the observed K~ are produced not too far away from the surface and therefore the 
multiplicity should not be too much influenced by the kaon properties at high density. Currently, 
one has approached this problem by using the same form of the mean field potential as for the 
with the parameters displayed in tabled The energy uj]^- [k — 0) decreases with increasing density 
and therefore one expects an increase of the K~ yield if the potential is switched on. This is indeed 
observed in the bottom figures of the figsllllandll2l The enhancement seen when the K~ potential 
is switched on is different in the different approaches. An important part of this difference comes 
from different extrapolations of the known K^N cross section, which diverges at small relative 
momenta to the ^/s < y/sthreshoU region. A solid theoretical basis for this extrapolation, which is 
necessary because the mass of the K~ decreases and hence the threshold is lowered, is not at hand 
presently. 

8. Lasting results 

One of the most challenging motivations for studying heavy ion collisions has been the search for 
the amount of energy which is needed to compress hadronic matter. This went under the slogan 
" Search of the nuclear equation of state" . Earlier ideas to use pions or the in- plane flow of baryons 
as an experimental signal for the compressibility of matter, have not led to unambiguous results. 
The pions turned out to be not sensitive because most of them are produced when the system 
has already expanded to a low density. The in-plane flow is small at densities reachable in these 
reactions and it depends on the pressure which is proportional to the density gradient. Nuclear 
surface properties, like density gradients, are in general much more difficult to simulate than bulk 
properties. 

It has been proposed to use for this purpose because at these subthreshold energies they 
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Figure 12. Final K~ transverse momentum distribution at b=l fm, \ycm < 0.5| and with 
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test the nuclear environment prior to production j42| . To get rid of many of the uncertainties it 
has been proposed by the KaoS coUaboration 0S| to use not the absolute yield but to compare a 
heavy with a light system where compression is virtually absent. Indeed, one has found that the 
excitation function of this ratio is rather sensitive to the compressibility of the nuclear potential 
and that the experimental results point strongly towards a nuclear potential with a compressibility 
K < 280MeV. Such a potential has usually been dubbed as the "soft equation of state". In the 
meantime, one has tested whether this result is robust, i.e. whether it remains valid if one varies 
the unknown or little known input such as the A lifetime or the AA'^ — > K'^ cross section with the 
conclusion that only a strongly density dependent production cross section may invalidate this 
result. Such a cross section, however, has never been proposed in elementary calculations. To which 
degree this result is stable against a change of the AN cross section is shown in the left panel 

of fig. 1131 where the different cross sections vary by a factor of 2. The results appear fairly robust. 
One should stress, though, that we find here that nuclear matter is relatively easy to compress. 
The compressibility is compatible with the results obtained by analyzing giant resonances, i.e. the 
compression of nuclei at very low temperature and around the equilibrium density. However, the 
reactions discussed here yield high densities and temperatures and the system contains mesons and 
resonances. Thus we test here a completely different region in the T-p plane. 



9. Conclusion 



Multiplicity, rapidity and transverse momentum distributions from different programs simulating 
heavy ion reactions at kinetic energies of around 1 AGcV have been compared. If the same inputs 
are used, the programs yield quite similar results but differ in details. Some of the differences 
are understood in terms of the implementation of different physical approaches in the programs, 
and others, less important, remain to be explained. Despite of these differences the programs 
agree on some important physical results: The pious are almost exclusively produced by resonance 
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Figure 13. Normalized ratio of the K'^ yield in Au+Au to C+C reactions as a function the 
beam energy for two different equations of states. Filled circles with errors represent the data of 
the KaoS collaboration 1431 . Other symbols with lines represent calculations of ref. 1211 on the 
left panel and of ref. 1191 on the right panel for either the AA^ — > A'+ cross section of ref. 1271 or 
for (o-(AA K+) = 3/4 a{NN K+)) employed earlier. 



decay. Subthreshold meson production is due to the same reaction as in free space and becomes 
possible because at least one of the collision partners has gained energy in collisions before. They 
programs furthermore agree on the stopping of matter, the dependence of the squeeze out 
on its in-medium potential and on the softness of the density dependence of the effective nuclear 
potential. 
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